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CN ■ 

(D ■ ABSTRACT 

PsJ I We improve upon the radiative, hydrostatic equihbrium models of passive circumstel- 

lar disks constructed by Chiang & Goldreich (1997). New features include (1) account 

^ ■ for a range of particle sizes, (2) employment of laboratory-based optical constants of 

2^ ■ representative grain materials, and (3) numerical solution of the equations of radiative 

■r^ij- , and hydrostatic equilibrium within the original 2-layer (disk surface + disk interior) 

^ , approximation. We systematically explore how the spectral energy distribution (SED) 

^^ I of a face-on disk depends on grain size distributions, disk geometries and surface den- 

^^ I sities, and stellar photospheric temperatures. Observed SEDs of 3 Herbig Ae and 2 T 

i-^ ' Tauri stars, including spectra from the Long Wavelength Spectrometer (LWS) aboard 

I ' the Infrared Space Observatory (ISO), are fitted with our models. Silicate emission 

5^ ■ bands from optically thin, superheated disk surface layers appear in nearly all systems. 

c/3 . Water ice emission bands appear in LWS spectra of 2 of the coolest stars. Infrared ex- 

. . , cesses in several sources are consistent with significant vertical settling of photospheric 

. ^^ I grains. While this work furnishes further evidence that passive reprocessing of starlight 

/\ ' by flared disks adequately explains the origin of infrared-to-millimeter wavelength ex- 

C^ ' cesses of young stars, we emphasize by explicit calculations how the SED alone does 
not provide sufficient information to constrain particle sizes and disk masses uniquely. 
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1. INTRODUCTION 

The energetics of the outermost regions of isolated disks surrounding T Tauri and Herbig Ae 
stars is dominated by passive reprocessing of central starlight. While many protostellar disks are 
actively accreting (see, e.g., the review by Calvet, Hartmann, & Strom 2000), effects of viscous 
dissipation on disk spectra manifest themselves most strongly in the immediate vicinities of the 
central stars, i.e., in the steepest portions of their gravitational potential wells. Simple scaling laws 
illustrate the relative importance of external irradiation vs. accretion luminosity. The local viscous 
luminosity per unit disk area decreases as 1/a^, where a is the stellocentric distance. By contrast, 
the flux of central stellar radiation striking the disk drops more slowly as (sina;)/a^, where a is the 
angle at which starlight grazes the disk surface. Vertical hydrostatic equilibrium normally ensures 
that disks flare outward such that a is a slowly increasing function of a for a S> i?* , where i?* is the 
stellar radius (see, e.g., Kenyon & Hartmann 1987). Hence, there is always a disk radius outside 
of which the energy from stellar illumination outweighs that of midplane accretion; in the extreme 
case that the central star derives its luminosity wholly from accretion, this transition radius is 
roughly 1 AU. The spectral energy distributions (SEDs) of young star /disk systems longward of 
~10^m should thus closely approximate those of passively heated disks, even when accretion is 
ongoing. 

Hydrostatic, radiative equilibrium models of passive T Tauri disks are derived by Chiang & 
Goldreich (1997, hereafter CG97). The passive disk divides naturally into two regions: a surface 
layer that contains dust grains directly exposed to central starlight, and a cooler interior that is 
encased and diffusively heated by the surface (Calvet et al. 1991; Malbet & Bertout 1991; CG97; 
D'Alessio et al. 1998). CG97 compute SEDs of passive disks viewed face-on and employ their model 
to satisfactorily fit the flatfish infrared excess and millimeter wavelength emission of the T Tauri 
star GM Aur. The optically thin, superheated surface layer is shown to be the natural seat of 
silicate emission lines (see also Calvet et al. 1992). 

In a second paper, Chiang &: Goldreich (1999, hereafter CG99) compute SEDs of passive T 
Tauri disks viewed at arbitrary inclinations. They point out that the spectrum of a nearly edge-on 
disk is that of a class I source for which i^Fi, rises from 2 to 10 ^m (Lada &: Wilking 1984; Lada 
1987). In general, class I sources are best described by a combination of an inclined, passively 
heated disk and a dusty bipolar outflow or partially evacuated envelope. The fraction of class I T 
Tauri spectra that represent limiting cases of simple, isolated, inclined disks is small (D'Alessio et 
al. 1999) but non-zero (CG99). 

This third paper in our series on passive protostellar disks extends our work in three directions: 



1. We refine equilibrium, 2-layer models of passive disks by (a) accounting for a range of par- 
ticle sizes, (b) employing laboratory-based optical constants of a suite of circumstellar grain 
materials, and (c) solving numerically the equations of radiative and hydrostatic equilibrium 
within our original 2-layer approximation. 

2. We systematically explore how the SED of a face-on disk depends on grain size distribu- 
tions, disk geometries and surface densities, and stellar photospheric temperatures. Physical 
explanations are provided for all observed behaviors of the SED. 

3. We employ our refined face-on models to fit observed SEDs of 3 Herbig Ae (HAe) and 2 
T Tauri stars. These observations include new spectra between 43 and 195 ;um from the 
Long Wavelength Spectrometer (LWS) aboard the Infrared Space Observatory (ISO) (Creech- 
Eakman, Chiang, van Dishoeck, & Blake 2000). The uniqueness of our fitted values of disk 
parameters is assessed, and evidence for emission lines from superheated silicates and ices is 
reviewed. 

The input parameters and basic equations governing our refined standard model are detailed 
in §2. Results, including a systematic exploration of how the SED varies in input parameter space, 
are presented in §3. Model fits to observations are supplied and critically examined in §4. There, 
we also compare our results to recent modelling efforts by Miroshnichenko et al. (1999). Finally, 
we summarize our findings in §5. 



2. Refined Model 

2.1. Input Parameters 

Table 1 lists the input parameters of our refined model. Figure 1 exhibits schematically the 
zones of varying grain composition in both the disk surface and disk interior. For ease of com- 
putation, and for want of hard observational constraints on the detailed compositions and optical 
properties of circumstellar grains, we limit ourselves to considering metallic iron (Fe, bulk density 
= 7.87 gcm^^), amorphous olivine (MgFeSi04, bulk density = 3.71 gcm^'^), and water ice (II2O, 
bulk density = Igcm^'^). These cosmically abundant materials span a wide range in condensation 
temperature (and therefore stellocentric distance), and in the cases of silicates and water ice, their 
existence is confirmed by spectroscopic observations (see, e.g., §4.3 of this paper). In reality, proto- 
stellar disks contain many more kinds of solid-state materials than we have incorporated. We have 
experimented with including additional grain compositions (e.g., graphite, organics, and troilite), 
but in no instance do we find our conclusions changed qualitatively. Our goal is not to be slavishly 
realistic but rather to highlight chief physical effects. 

Thus, where local dust temperatures (= gas temperatures) fall below T*^^ ~ 150 K, the grains 
are taken to be spheres of amorphous olivine mantled with water ice. For simplicity, the thickness 



Table 1. Input Parameters of Refined Model 



Symbol 



Meaning 



Standard Value 



So 
P 

tto 

H/h 
Qi 



rpzron 

sub 

rpsil 

sub 

'-pice 

sub 

Mdisk' 



Stellar Mass 

Stellar Radius 

Stellar Effective Temperature 

Surface Density at 1 AU 

— d log S/d log a 
Outer Disk Radius 

Visible Photospheric Height / Gas Scale Height 

— d log N/d log r in Interior 

— d log N/d log r in Surface 
Maximum Grain Radius in Interior 
Maximum Grain Radius in Surface 
Iron Sublimation Temperature 
Silicate Sublimation Temperature 
H2O Ice Sublimation Temperature 
Total Disk Mass (Gas + Dust) 



O.SM© 

2.fyRQ 

4000 K 

10^ gcm^^ 

1.5 

8600i?* = 100 AU 

4.0 

3.5 

3.5 

1000 ;um 

1/im 

2000 K 

1500 K 

150 K 

0.014 Mq 



^The total disk mass is not an explicitly inputted parameter but is derived from 
So, p, Co, and the inner disk cut-off radius, Oj = 2i?*. 



of the water ice mantle, Ar, relative to the radius of the olivine core, r, is held constant. Where 

pice „„J rpsil 
-sub ^^'^ ^sub 



local temperatures fall between T**;? and T**', ^ 1500 K, only the pure olivine cores are assumed to 



remain. In innermost disk regions where local temperatures fall between T^^'^ and Tj^^"^ ~ 2000 K, 
the grains are taken to be spheres of metallic iron. 

The iron or silicate cores in the disk surface (interior) possess a power-law distribution of radii 
r Detween t^yi'Iyi ano. T^yi^x^s y'^max^i) • 

dN oc r~«-(') dr , (1) 

where dN is the number density of grains having radii between r and r + dr. Variables subscripted 
with "s" denote quantities evaluated in the disk surface, while those subscripted with "i" denote 
quantities evaluated in the disk interior. Our standard value of Qi = Qs = 3.5 places most of the 
geometric surface area in the smallest grains and most of the mass in the largest grains. In practice, 
Vmin is fixed at 10"^ /^m, while rmax,i, rmax,s, Qi, and Qs are free to vary. Generally rmax,s < rmax,i 
since large grains tend to settle quickly out of tenuous surface layers (see §3.3 of CG97). All of the 
cosmically abundant iron is assumed to be locked within grains. Following Pollack et al. (1994), 
we take 50% of the cosmically abundant oxygen to be locked in H2O ice. Values for all cosmic 
abundances are obtained from Allen (2000), except for the abundance of oxygen which is taken 
from Meyer, Jura, & Cardelli (1998). Together, these assumptions yield a fractional thickness, 
Ar/r, for the water ice mantle equal to 0.4. 

Optical constants for amorphous olivine are obtained from the University of Jena Database 
[http://www.astro.uni-jena.de; see also Jager et al. (1994)]. Longward of 500 /im where optical data 
for silicates are not available, the complex refractive index (n + ik) for glassy olivine is extrapolated 
such that n(A > 500 ^um) = n(500;um) and k{X > 500 ^um) = k{500fim) (A/500 /im)"^ Optical 
constants for pure crystalline H2O ice are taken from the NASA ftp site (ftp : climate.gsfc.nasa.gov 
/ pub / wiscombe / Refrac_Index / ICE /; see also Warren 1984). Though employing the constants 
for a cosmic mixture of amorphous ices (H2O : CII3OII : CO : NH3; see Hudgins et al. 1993) 
would be more appropriate, we nonetheless adopt the data for pure II2O ice because the latter 
are available over all wavelengths of interest, from the ultraviolet to the radio, whereas the former 
are not. One consequence of using the constants for crystalline (213-272 K) water ice as opposed 
to amorphous (~100 K) ice is that spectral features due to translational lattice modes at 45 and 
62;um are slightly underestimated in width and overestimated in amplitude (see, e.g., Hudgins et 
al. 1993). Optical constants for metallic Fe are obtained from Pollack et al. (1994). 

The inner cutoff radius of the disk, Oj, is fixed at 2i?*. For T Tauri stellar parameters, this 
radius coincides with the distance at which iron grains in the surface layer attain their sublimation 
temperature. For the hotter HAe stars, the iron condensation boundary occurs at a ~ 14-30 -R^,. 
Inside the iron condensation radius, the disk may still be optically thick to stellar radiation even if 
dust is absent. Opacity sources include pressure-broadened molecular lines and Rayleigh scattering 
off hydrogen atoms [see the appendix of Bell & Lin (1994)]. For simplicity, when modelling HAe 
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Fig. 1. — Schematic of zones of grain compositions (iron -^ amorphous olivine -^ amorphous olivine 
mantled with water ice) for both the disk surface and interior. The numerous vertical lines drawn 
in the surface indicate that there we account for how grains of different sizes sublimate at different 



stars, we employ a 1-layer blackbody disk that extends from the iron condensation boundary to 



2.2. Grain Absorption Efficiencies and Opacities 

The grain emissivity, e(r, A), is equal to its absorption efficiency and is calculated using Mie- 
Giittler theory (Bohren & Huffman 1983; see their subroutine BHCOAT.F).^ Figure 2 displays 
absorption efficiencies for our ice-silicate, silicate, and iron spheres having three representative 
sizes. The emissivity index in the Rayleigh limit, (3 = dine/ din i/, equals 1.64, 2.00, and 0.50 for 
the three compositions, respectively. 

Well-known resonances at A < 20 /um include the 0-H stretching (3.1, 4.5 ^m) and H-O-H 
bending (6.1 /im) modes in water, and the Si-0 stretching (lO^um) and 0-Si-O bending (IS^um) 
modes in silicates. 

In crystalline water ice (and crystalline silicates), optically active modes of vibration longward 
of '-^10/um are "intermolecular translational" or "intermolecular rotational." These involve collec- 
tive movement of a molecule or a unit cell with respect to other molecules/unit cells in the lattice. 
The strengths, positions, and widths of these modes are more sensitive to the presence of chemical 
impurities and to long range order in the solid (i.e., its degree of crystallinity, or "allotropic state") 
than those of fundamental stretching and bending modes at shorter wavelengths. In principle, 
these intermolecular modes provide information on the annealing history of initially amorphous, 
interstellar material in the relatively high density and high temperature environments of circum- 
stellar disks. The intermolecular translational modes in water ice evident in Figure 2 are located 
at 45, 62, and 154 ^m (Bohren & Huffman 1983, see page 278; Bertie, Labbe, & Whalley 1969, see 
their Figures 4 and 11). All of these features are positioned within the wavelength range of the 
Long Wavelength Spectrometer aboard ISO. Note also the blending of the 12 /xm intermolecular 
rotational band in water ice (Bohren & Huffman 1983) with the Si-0 silicate stretching mode at 
10 /im. 

Oscillatory behavior in Figure 2 near the onset of the Rayleigh limit (27rr/A ~ 1) reflects so- 
called "ripple structure" that arises from our use of perfectly spherical particles (Bohren & Huffman 
1983); we expect real- world deviations from sphericity to smooth out this artificial behavior. 

In the disk interior, the opacity is given by 

^.(A) = -/ —r^e{X^r)dT, (2) 

where pt is the total density of gas and dust. Figure 3 displays Kj for our distribution of ice-silicate 



®Our grain emissivity, e, equals Qaba in the notation of Bohren & HufTnian (1983). 
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Fig. 2. — Emissivities of ice-silicate (H2O / amorphous olivine), silicate (amorphous olivine), and 
iron grains having three representative core sizes. The thickness of the water ice mantle relative 
to the radius of the olivine core is Ar/r = 0.4. Resonant features include the 0-H stretching (3.1, 
4.5 //m) and H-O-H bending (6.1 /xm) modes in water ice; the Si-0 stretching (10 /um) and 0-Si- 
O bending (18 /im) modes in silicates; and the intermolecular translational (45, 62, and 154 ^um) 
modes in water ice. Oscillatory behavior near the onset of the Rayleigh limit (27rr/A ~ 1) reflects 
"ripple structure" arising from our use of perfectly spherical particles. 



and silicate particles in cosmic abundance gas. 



2.3. Basic Equations 

In the surface layer, dust grains are directly exposed to stellar radiation. The grains attain an 
equilibrium temperature 

where {e{r))T is the emissivity of a grain of radius r averaged over the Planck function at temper- 
ature T, and ~ 1/2 is the fraction of the stellar hemisphere that is seen by the grain. 

In our two-layer formalism, exactly half of the radiation reprocessed by the surface layer escapes 
directly into space. The remaining half is directed towards the disk interior. Of this latter half, a 
fraction, 1 — e~^''^»'«, is absorbed by the disk interior, where S is the disk surface density and (kj)s 
is the opacity of the disk interior averaged over the spectrum of radiation from the surface. ^'^ In 
radiative balance, 

^^l_^-^{^.h^ f^\ Ttsma = {l-e-''^^^^^)T,\ (4) 

where Tj is the temperature of the disk interior and («;«)« is the interior opacity averaged over the 
spectrum of radiation from the disk interior. ^^ Note that Tj is the common temperature of interior 
grains of all sizes which are assumed to have thermally equilibrated with one another. Here a is 
the angle at which stellar radiation strikes the surface: 

^ (d\nHH\ ^ H / 4 RA 

a ~ arctan I — — arctan h arcsm (5j 

\ dlna a J a ySir a J 

[cf. equation (5) of CG97]. The height of the disk photosphere, H, is assumed to be proportional to 
the vertical gas scale height, h, with a fixed constant of proportionality equal to 4 for our standard 
model. In reality, when dust and gas are well-mixed in interstellar proportions, the ratio H/h 
decreases slowly from ~5 at 1 AU to ~4 at 100 AU. In modelling observed SEDs, we will allow 
H/h to be a fitted constant parameter. In hydrostatic equilibrium. 



^"in practice, we perform this average over a Planck function evaluated at the temperature of the most luminous 
grains in the surface. For our assumed size distribution, these dominant grains typically have radii r ~ 0.5 /im; surface 
grains having r « 0.2-1 /im are responsible for absorbing ~50% of the incident stellar radiation. 

^^We perform this average over a Planck function evaluated at temperature Tj. 
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Fig. 3. — Mass absorption coefficients for our standard model distributions of ice-silicate and silicate 
particles in solar abundance gas. The same resonant features found in Figure 2 are seen here. We 
have smoothed these curves to suppress so-called "ripple structure" at A > 1000 //m that arises 
from the sphericity of our particles having r ~ 1000 fiui. Note that these curves are sensitive to our 
chosen qi and rmax,i- 
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(6) 



where Tc = GM^^Hg/kR^ and ^u^ = 3 x 10 ^'^ g 



Equations (4) and (6) are two equations for the two unknown functions, H{a) and Ti{a). 
Substitution of (6) into (4) yields an algebraic equation for Tj and the slowly-varying flaring index, 
7 = din H/ din a. We solve this equation for Tj(a) and 7(a) numerically on a logarithmic grid in 
stellocentric distance. Our procedure is described in detail in the Appendix. 

The SED of the disk equals the sum of emission from the disk interior, 

4TTd'^XFx,i = Svr^A / " Bx{Ti) (1 - e-^''')ada, (7) 

J ai 

and from the surface layers (above and below the disk midplane). 



rao 

47rd2AFA,, = 87r2A(l + e-^^0 / Sx{l - e-^')ada , 

J ai 



(8) 



where d is the distance to the source. The source function in the surface, S\^ is the Planck function 
averaged over the ensemble of effective grain cross-sections in the surface layer: 



2j;-^-Bx{Tas)^r^e{\T)dr 

■^ ' rnvn '^ 

^ ' mm ai 






In equation (9), the factor of 2 is inserted so that exactly half of the incident radiation is reprocessed 
by the surface layer. The normal optical depth of the surface layer, tx, is similarly averaged: 



3. Results 

3.1. Flaring Index 

Figure 4 displays the behavior of the flaring index, 7 = dlnH/dlna. We conceptually divide 
the disk into three annular regions, as was done in CG97 (see their §2.3.2). In the region marked "I," 
the disk interior is opaque to both its own reprocessed radiation and to radiation from the surface. 
Here 7 increases from its flat disk value of 1.125 ~ 9/8 to its asymptotic value of 1.275 ~ 9/7 as 
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the first two terms on the right-hand side of equation (5) gradually dominate the last term. In 
region "II," the disk interior remains opaque to radiation from the surface, but is optically thin 
to its own reprocessed radiation. Here 7 steeply rises with a because grains in the disk interior 
equilibrate at relatively high temperatures to compensate for the relative inefficiency with which 
they re-radiate the incident energy. Finally, in region "III," the interior is transparent to radiation 
from the surface (i.e., S(«;j)s < 1); the inability of the interior to absorb the incident energy causes 
7 to decrease. 



3.2. Disk Temperatures 

Figure 5 exhibits temperature profiles for the surface (Tds) and for the interior (Tj). The 
temperatures of grains in the surface layer vary slightly with their sizes; in Figure 5, we have 
chosen to plot T^s for the size bin containing the most luminous grains, i.e., the logarithmic size 
interval that absorbs the greatest fraction of incident stellar radiation. For our choices of grain 
composition and size distribution, these dominantly absorbing grains have radii r ~ 0.1-0.7 /um. 
For reference, we also overlay in Figure 5 the temperature of an imaginary blackbody sphere, Tbb, 
which is naked before half of the stellar hemisphere. 

These temperature profiles are largely similar to those found in the simpler model presented 
in CG97, and help to justify the approximations made there. At a given distance, the surface is 
hotter than the interior by a factor of ~3. Though the interior in region "II" (see §3.1 above) is not 
radially isothermal as was found in the cruder analysis of CG97, Ti{a) does flatten slightly at these 
distances, as expected. Deviations from a single power-law behavior for T^sia) arise from structure 
in e{X). For example, T^s declines slightly more steeply with a for 0.5 ^ oau ^ 2 because grains 
at these distances cool relatively efficiently through silicate resonances at 10-20 //m. 



3.3. Refined Standard SED 

The SED for our refined standard model is displayed in Figure 6. The result shown here and 
that of CG97 (see their Figure 6) differ primarily in the emission from the superheated surface. 
Here we have accounted in detail for the optical properties of a few likely circumstellar grain 
materials. Solid-state resonances of superheated dust grains residing in disk surface layers appear 
in emission. These include vibrational modes in silicates at 10 and 18 ^m, and lattice translational 
modes in water ice at 45 and 62 ^m. Emission lines from vibrational resonances in ice shortward of 
10 /im are absent; ice is not present in the disk surface inside 6 AU.^^ Note the dearth of emission 



^^It is conceivable that vibrational resonances in water shortward of 10 /im may still appear in disk spectra if silicate 
particles are hydrated. Such signatures have been observed in spectra of CI chondritic meteorites (McSween, Sears, 
& Dodd 1988). 
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Fig. 4. — The flaring index, 7 = din H/dln a, for our refined standard model. As was done in CG97 
(see their §2.3.2), we divide the disk into three annular regions depending on the optical depth of the 
disk interior. In region "I," the interior behaves as a blackbody; 7 increases from its flat disk value of 
1.125 ~ 9/8 to its asymptotic flared value of 1.275 ~ 9/7 as the disk thickness becomes increasingly 
larger than the stellar radius. In region "II," the disk interior becomes optically thin to its own 
reprocessed radiation; 7 increases as interior grains enhance their temperatures to compensate for 
the inefficiency with which they re-radiate. In region "III," the interior is transparent to radiation 
from the surface and cools quickly with increasing distance, causing 7 to decrease. 
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Fig. 5. — Temperature profiles for the surface and for the interior in our refined standard model. 
The temperatures of grains in the surface layer depend on their sizes; here, the curve marked Tds 
represents the size bin, r ~ 0.5 /um, containing the most luminous grains. The discontinuity in Tds 
at a ~ 6 AU marks the water condensation boundary in the surface, outside of which H2O ice coats 
silicate cores; the discontinuity in T^s at a ~ 0.05 AU marks the silicate condensation boundary in 
the surface. For reference, the temperature of a spherical blackbody which is naked before half of 
the stellar hemisphere is shown as a dashed line. 
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Fig. 6. — Spectral energy distribution for our refined standard model. Resonant bands of super- 
heated dust grains in disk surface layers appear in emission. These include vibrational modes in 
silicates at 10 and 18 /im, and lattice translational modes in water ice at 45 and 62 /xm. Absent are 
emission lines from vibrational resonances in ice shortward of 10 /um; ice is not present in the disk 
surface inside 6 AU.^^ Note the dearth of emission from the surface between 2 and 8;um; amor- 
phous olivine (MgFeSi04) grains having r < 1/im are relatively transparent at these wavelengths 
(see Figure 2). 
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from the surface between 2 and 8/xni; amorphous ohvine (MgFeSi04) grains having r < 1 iiui are 
relatively transparent at these wavelengths (see Figure 2). The signature of this "silicate transparent 
region," however, is largely masked by emission from the optically thick interior at small radius. 
Furthermore, our computational experiments demonstrate that including other types of particles 
in the surface layer such as troilite (FeS) serves to further fill in this transparent region. The broad 
peak in disk surface emission near A ~ 1.5 ^m arises from our pure iron particles and iron impurities 
in our olivine particles. 



3.4. Dependence of SED on Input Parameters 

In Figures 7 and 8, we explore the dependence of the SED on our input parameters. In each 
panel, we vary the indicated parameter(s) while fixing all other parameters at their standard model 
values. All the variations are easily understood. We observe the following behavior: 

1. Millimeter- wave fluxes are most sensitive to So, p = —dlnT,/dlna, rmax,i, sand qi = —dlnN/dlnr. 
The former two variables determine the amount of mass in the cool disk interior at large ra- 
dius; the latter two variables affect the millimeter-wave opacity in the disk interior. 

2. Millimeter-wave SEDs for rmax,i = 1 and 10 //m are identical; these two cases belong to the 
Rayleigh limit, in which absorptive cross-sections are proportional to grain volume (Bohren 
&: Huffman 1983). In this limit, millimeter- wave opacities are independent of how the total 
condensable mass in water and silicates is distributed with particle size. 

3. As rmax,i increases from 10 to 100 /um, the opacity at A = 100-600 ^m in the disk interior also 
increases, thereby enhancing emission from the disk interior at those wavelengths. A similar 
effect is seen as r^ax,! increases from 100 to 1000 //m. 

4. For all values of rmax,i considered, the spectral index of the SED at A = 2-4 mm equals 
77-2-4 = dln{iy'Fi^)/dln u = ?> + f3es = 4.6. The value of f3cs = 1-6 equals the value of fi for our 
ice-silicate grains, indicating that radiation at these wavelengths emerges from optically thin 
material. 

5. For qs = qi = 2.5, most of the geometric cross-section in our dust size distribution is concen- 
trated in the largest grains {rmax,s ~ 1 Mni in the surface and rmax,i ~ 1000 ^m in the interior). 
Compared to the standard model, silicate emission features from the surface at mid-infrared 
wavelengths are weaker because the Rayleigh limit does not apply for these large grains. At 
millimeter wavelengths, n2-4 = 4.3 =^ /3eff = 1.3 < (3ice~sa = 1-6; a substantial fraction of the 
mm-wavelength emission arises from the disk interior made optically thick by the increased 
number of mm-sized particles. Note also the relative dearth of emission from A ~ 80-400 /um 
compared to the standard model, caused by fewer numbers of r ~ 15-70 ^m sized particles 
in the disk interior and a concomitant loss of interior opacity at these wavelengths. 
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Fig. 7. — Dependence of SED on input parameters So, p, rmax,i, and (/j = Qs- In each panel, the 
indicated parameter takes on several values while all other parameters are fixed at their standard 
model values. Grain size indices qi and qs are varied simultaneously for compactness of presentation. 



18 



(7J 


C\2 


T— 1 


o 


O) 


CO 


C\2 


T— 1 


o 


O) 


CO 


CO 


CO 


CO 


C\2 


CO 


CO 


CO 


CO 


cv 



:l 



< 



T—\ 


_yr 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


—r 


~y^ 


r 1 1 1 1 


1 1 1 1 II 


1 1 1 1 1 1 1 1 




:\ ^ c\3 -H 


- 


-\ 


. 


D 


D D - 








- 


\ 


< 


< < -_ 


o 


= ^^5- 


1 ; 
1 I 


- 


E 


y 


o 


o o - 


T— 1 




1 : 

1 \ 






^ 


o 


o o " 






1 : 





_ 


\ 


T— 1 


C\2 CO J 


o 




- 


^ 








■ 1 — 


o 


~ \\ ■. 


_ 


=- 


\ 
\ 






; 1 


T— 1 


\x'-. 






%. 






: 1 




\x"-. 




~ 


\i^ 










v^._ 





_ 


"• 


V 




: 1 












^^ 




o 


^^^ 


- 


- 




^^*-^ 


- 


o 


E ,rl ^'^'^^^^ 


- 


= 


o 


^*^=. 


"*»=>>. 


o 


~ ^\ ^^'^'"'^i^ 


- 




cd 




^'^^'^^'^^-^^ - 


T— 1 


: X 


■*&*. 


: 






^^ 




=1 1 1 1 1 1 1 1 1 1 1 1 1 1 III 


,i. 


= 






1 l" 




ll\l \l r-l 1 l>l^ 1 1 1 1 1 1 III 111 


1 1 


-y 


' ' ' ' 


III II 


1 III II 


T 1 


E ^ \ ■• ( 


- 






s 


e - 






— 


- 


\ 


6 i - 




= \ \ '••. \ o o o o 
~^^) <■■■'■ ^ oooo 


- 


- 


\ 




a.^ - 


o 


~ 


=_ 


^ 


o 




T— 1 


>- ;-■■ > oooo 

: "'A \ \ ^ CD 00 ^ 


" 




y 


T— 1 


— 1 d - 




— 


- 


\ 






— 


o 
o 

T 1 




; 1 1 
: 1 
: 1 


- 


=- 


\ 


1 

1 
1 
1 




; _ 




xX 


" 


— 


N^ 




~ 




: ^^ 


— 


- 


m 




— 


o 


"^ 


- 


^ 


- 




- 


o 


= ^^ 


- 


^ 


X 




■v^ 


o 




- 




CO 

s 

1 1 1 1 




^"^^^^ - 


T— 1 


: ^ 

=1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 


=1 


1 1 1 1 1 1 


1 1 1 1 1 1 1 1 



^ 



CO 



C\2 



ISl 

T 1 ■^^— ^ 

O 

^ o 



CO 



C\2 



LO^coc\2^oa:) cooj^oo) 

C0C0C0C0C0C0C\2 C0C0COC0C\2 

(s/§ja) ^j/i^pixt' DOT 

Fig. 8. — Dependence of SED on input parameters T*, H/h, rmax,s, and Oq- In each panel, the 
indicated parameter takes on several values while all other parameters are fixed at their standard 
model values. 



19 



6. For Qs = Qi = 4.5, most of the condensable mass is concentrated in the smahest grains 
(rmin ~ 0.01 ;um). Mihimeter-wave opacities and therefore fluxes are lower than for our 
standard model. Emission at A < 100 ^m does not differ from standard model results because 
the Rayleigh limit still applies for the most luminous grains in the surface; shapes of the 
emissivity curves for silicate and ice-silicate grains are independent of r in the Rayleigh limit 
(see Figure 2). 

7. The radial locations of condensation boundaries in disk surface layers move outward approx- 
imately as T^ ~ T^. Consequently, as T^ increases, surface emission from water ice 
diminishes noticeably. 

8. Reducing the height of the disk photosphere by reducing the scaling parameter H/h (thereby 
crudely modelling the effects of vertical settling of dust) lowers the amount of stellar radiation 
intercepted and reprocessed by the disk. Emission at A < 200 /im scales nearly linearly with 
H/h. Radiation at these wavelengths arises from the optically thick interior (region "I") and 
from the optically thin surface; for both regimes, vFi, scales as sin a which scales approxi- 
mately as H/h. Radiation at A > 200 ;um — the Rayleigh-Jeans regime — is less sensitive to 
H/h; here the approximate scaling relation reads uF^^ oc Tj oc (sina)*^'^^ oc (H/h)^''^^. 

9. Surface SEDs vary negligibly with rmax,s- For our standard slope of the size distribution 
(qs = 3.5), silicate and ice-silicate particles having radii r* ^ 0.5 ^m absorb the bulk of 
the radiation from the T Tauri star; surface grains having r ~ 0.2-1 fim are responsible for 
absorbing '^50% of the incident stellar radiation. If rmax,s > ^*, the SED is unchanged 
because those grains having r ^ r^, are insufficiently numerous to be significant absorbers of 
radiation. If r„iax,s < f*, the surface SED remains unaltered because the Rayleigh limit still 
obtains. 



4. Fitting SEDs of T Tauri and HAe Stars 

Our sample comprises 3 HAe and 2 T Tauri stars that are (1) not known to harbor stellar 
companions, and (2) not known to drive jets or to be surrounded by massive, > 500AU-scale 
nebulosities that are better described by spherical envelopes rather than by fiattened disks. In 
order of decreasing stellar effective temperature, the sample stars are MWC 480 (= HD 31648), 
HD 36112 (= MWC 758), CQ Tau (= HD 36910), LkCa 15, and AA Tau. For 4 of our sources 
(MWC 480, HD 36112, CQ Tau, and AA Tau), medium-resolution (AA = 0.2 ^m) ISO LWS 
spectra between 43 and 195 ^m are available. More detailed descriptions of the ISO data, including 
how they were reduced and what gas phase spectroscopic information they contain, are presented 
elsewhere (Creech-Eakman, Chiang, van Dishoeck, & Blake 2000). All of our sources were too weak 
to be observed by the ISO Short Wavelength Spectrometer. 

The code that computes our refined standard model is restricted to calculating SEDs for face- 
on disks. A detailed study of how the SED varies with inclination, i, has been given previously 
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(CG99). As described there, non-zero inclinations affect the infrared SED in 2 principal ways: at 
moderate i, by introducing a cosine i variation in emission arising from the optically thick interior, 
and at extreme i, by blocking radiation from the disk at small radius via the intervening flared disk 
at large radius. The second effect on the infrared SED is negligible for our sample stars. Visual 
extinctions range from Ay = 0.3 mag (MWC 480) to 2 mag (CQ Tau); these modest values imply 
that inner disk regions are not significantly occulted by flared outer disk edges at mid-to-far infrared 
wavelengths. We estimate that the first effect introduces at most a factor of 2.5 overestimation in 
our computed fluxes between 2 and 8 ^m where the SED is dominated by emission from the optically 
thick interior. For example, the disk inclination for CQ Tau has been independently estimated to be 
~66°, based on photometric and polarimetric variability at visible wavelengths (Natta & Whitney 
2000). Such "UXOR- type" phenomena has been interpreted to arise from clumps of dust in the 
flared disk surface sporadically obscuring our line-of-sight to the central star. 

For each source, a model SED is fitted to the ISO LWS scan (if available), millimeter wavelength 
fluxes, and ~3-25 /um photometric data. In 3 of the sources (AA Tau, CQ Tau, and MWC 480), 
ISO fluxes are greater than corresponding IRAS (Infrared Astronomical Satellite) fluxes at 60 
and 100 ^m by factors of ~2-3. The origin of the discrepancies is not known. Where there are 
discrepancies, preference is given to the ISO LWS data for which the beam area is ~2 times smaller 
than that of IRAS. For HD 36112, there is excellent agreement between ISO and IRAS. Preference 
is given also to the central portions of the ISO scans between 50 and 170 ^m where individual 
detectors overlap in wavelength coverage and measured fluxes are consequently more reliable. 

In fitting the SEDs, we fix Qs = 3.5. Smaller values {qs < 3) seem unlikely since they would 
imply that the largest grains, which tend to settle out of surface layers most quickly, dominate the 
geometric cross-section. The SED is insensitive to larger values {qs > 4), as shown in §3.4. We also 
fix rmax,s = l^m and p = 1.5 for all models. Section 3.4 demonstrates that the SED is insensitive 
to rmax,s once Qs is fixed at 3.5, and that Sq and p affect the SED in similar ways. 

Observations and fitted theoretical models are displayed in Figures 9-13. Table 2 contains the 
fitted parameters for our sample. The fits are intended to be illustrative; no attempt is made to 
minimize fit deviations in a formal, statistical sense. The results of such an analysis would not 
be very meaningful anyway, since the SED tends to be degenerate with respect to simultaneous 
changes in several of the parameters, as we discuss in §4.1. In any case, our fitted outer disk 
truncation radii for CQ Tau, HD 36112 (MWC 758), and MWC 480 are in accord with upper limits 
based on A = 2.7mm continuum images taken by Mannings Sz Sargent (1997). 

With few exceptions, the agreement between models and observations is good to within a 
factor of 2, and serves as further evidence that simple reprocessing of central starlight by flared 
disks adequately explains the infrared excesses of these systems. We are aware that this is not 
the conclusion of Miroshnichenko et al. (1999, hereafter MIVE99), who require the presence of 
tenuous envelopes having radii ~1000AU to heat embedded HAe disks to temperatures greater 
than those predicted by the classical T oc a^^'^ law. This idea was first elucidated by Natta (1993) 
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Fig. 9. — Refined 2-layer model fitted to data for MWC 480. Dotted lines denote the contribution 
from surface layers. Dashed lines denote the contribution from the disk interior. The inset plot 
with linear ordinate magnifies the spectral region spanned by ISO LWS. Photometric data are taken 
from Mannings & Sargent (1997), the color-corrected IRAS Point Source Catalog (1988), Sitko et 
al. (1999), Thi et al. (2000), and Qi (2000). 
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Fig. 10. — Refined 2-layer model fitted to data for HD 36112. Photometric data are taken from 
Mannings & Sargent (1997), tlie color-corrected IRAS Point Source Catalog (1988), and Thi et 
al. (2000). 



23 



1000 



A (yu-m) 
100 10 



OT 






w 

o 
o 



35 
34 
33 
32 
31 
30 
29 
28 



1 1 1 1 1 1 I I — I — I 1 1 1 1 1 I I — I — I 1 1 1 1 1 I I — I — I 1 1 1 1 1 I I 

CQ Tau 

'^max.i - 1000 Mrn q, = 3.5 

Eq = 2000 g/cm^ a^ = 180 AU 

H/h = 5.0 



\ 

\ 

\ 
- \ 
\ \ 




'l 1 1 1 1 1 1 1 1 1 1 I I I I I I I I I I I 1 I I I — I — I — I — I — I — I — I 1 — "n r 

Linear ordinate 
(units of Lq) 




III I I I I I I I I I I I I I L 



2 

1.5 
1 
0.5 



- 



200 150 



100 

A (yLim) 



50 



Fig. 11. — Refined 2-layer model fitted to data for CQ Tau. Photometric data are taken from 
Mannings & Sargent (1997), the color-corrected IRAS Point Source Catalog (1988), and Thi et 
al. (2000). 
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Fig. 12. — Refined 2-layer model fitted to data for LkCa 15. Photometric data are taken from the 
color-corrected IRAS Point Source Catalog (1988), Thi et al. (2000), and Qi (2000). 
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Fig. 13. — Refined 2-layer model fitted to data for AA Tau. Photometric data are taken from 
Beckwith et al. (1990), Beckwith & Sargent (1991), and Dutrey et al. (1996). 
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Table 2. Fitted Parameters^ of Herbig Ae and T Tauri Star/Disk Systems 



Parameter MWC 480 HD 36112 CQ Tau LkCa 15 AA Tau 



T, (K) 


8890 


8465 


7130 


4395 


4000 


R*{Rq) 


2.1 


2.1 


1.27 


1.64 


2.1 


M,{Mq) 


2.3 


2.2 


1.7 


1.0 


0.67 


d (pc) 


140 


150 


100 


140 


140 


So(gcm-2)'= 


8000 


1000 


2000 


9000 


1500 


ao (AU)^ 


100 


250 


180 


200 


250 


H/h 


1.7 


1.5 


5.0 


1.0 


3.8 


q^'' 


2.8 


3.5 


3.5 


2.5 


3.5 


rmax,i il^T^T 


1000 


1000 


1000 


3000 


1000 


MDISK(M0)d 


0.11 


0.02 


0.04 


0.18 


0.03 



H{ao)/ao '^ 0.13 0.16 0.45*' 0.09 0.58 

^For all sources, we fix qs — 3.5, rmax,s = 1 /^ni, and p — 1.5. See text for rationale. 

Stellar parameters and distances for HAe stars are taken from Mannings & Sargent (1997), except for 
CQ Tau for which R, and d are normalized to the Hipparcos distance. Stellar parameters and distances for 
T Tauri stars are taken from Beckwith et al. (1991) and Webb et al. (1999). 

"^The continuum SED is largely degenerate with respect to simultaneous changes in So, rmax,i, <li, and 
Oq. The values shown here are not uniquely constrained. 

Total masses (in gas and dust) and maximum aspect ratios of fitted disks are derived quantities and not 
input parameters. As discussed in §4.1, the total disk mass depends sensitively on the a priori unknown 
millimeter- wave opacity and could vary by an order of magnitude. 

''Our fitted maximum aspect ratio for CQ Tau and the modest visual extinction to this star (2 mag) imply 
a disk inclination i < 66° (i — for face-on views). This is consistent with the independent estimation of 
i ^ 66° by Natta & Whitney (2000) based on UXOR-type variability at visible wavelengths. 
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in the context of flat spectrum T Tauri stars. While some HAe systems do exhibit large scale 
nebulosities whose sizes as functions of wavelength vary in accord with the calculations of MIVE99, 
we disagree with the statement that disks are incapable of explaining the SEDs of HAe systems 
without heating from an envelope. Non-classical temperature laws follow naturally from passive 
reprocessing of starlight by hydrostatically flared disks (e.g., Kenyon & Hartmann 1987; CG97; this 
paper). 

Moreover, MIVE99 claim that radiation from disk surface layers as envisioned by CG97 do 
not contribute substantially to SEDs of HAe stars MWC 480 and CQ Tau. However, in Figures 9 
and 11, our model fits to these 2 sources indicate that (1) emission from optically thin disk surface 
layers dominates emission from the optically thick disk interior between 10 and '^50;um, and (2) 
surface layer emission naturally explains the observed infrared excesses, in particular the presence 
of silicate and ice emission bands (see §4.3 for a more complete discussion of these bands). 

Certainly the SED alone does not furnish sufficient information to uniquely constrain the 
geometry of dust surrounding HAe stars. We agree with MIVE99 that the presence of envelopes is 
an important possibility to consider for all HAe stars, and that for sources such as AB Aur their 
existence is persuasively implicated by imaging data at multiple wavelengths (see §2 of MIVE99). 
We simply wish to emphasize here that the SED alone does not rule out isolated disks heated by 
their central stars, and that envelopes are not the only way to achieve extra heating of the disk. 



4.1. Degeneracy between Disk Mass and Grain Size Distribution 

As might be gleaned from Figures 7 and 8, the values for Sq, p, rmax,i, Qi, and Oo presented in 
Table 2 cannot be uniquely constrained by the continuum SED alone. We display one degenerate 
combination in Figure 14, where 2 models using 2 different sets of parameters are fitted to the 
observed data for HD 36112. We feel that the fits are of comparable quality, given the crudeness 
of our 2-layer model. It follows that the total disk mass (in dust) is highly uncertain; in fact, 
our two models for HD 36112 differ in their total dust mass by a factor of ~6. One relies on the 
optical thinness of the disk interior at millimeter wavelengths to measure the disk mass in dust: 
-^DISK ex -Fmm oc Tmrn oc Skj. The problem is that the SED alone cannot disentangle S from Kj, 
which in turn depends on qi and rmax,i- 

The two models in Figure 14 yield significantly different results for the mass distribution in 
space, however. Upcoming spatially resolved observations at millimeter wavelengths will help to 
break this degeneracy between grain size (i.e., millimeter- wave opacity) and disk surface density 
(Beckwith, Henning, & Nakagawa 2000). 
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Fig. 14. — Degeneracy between surface density and grain size for HD 36112. Two models having 
different sets of parameters are fitted to tlie same dataset. Tfie disk mass, -MdisKj is not an explicitly 
inputted parameter but is derived from Sq, "p = 1.5, Co, and a-i = 2R^. Smaller grain sizes (i.e., 
smaller millimeter-wave opacities) may be traded for larger disk masses to achieve similar mm-wave 
optical depths and fits of comparable quality. Note, however, that both models require H/h = 1.5 
to match the overall level of infrared excess at A < 100 lira.. 
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4.2. Evidence for Dust Settling 

The one disk parameter which appears to be most uniquely constrained is H/h, the height of 
the disk photosphere in units of the gas scale height. Its value is roughly proportional to the overall 
level of infrared excess at A < 100 jiva.. CQ Tau exhibits H/h = 5.0, a value appropriate for gas and 
dust that are well-mixed in interstellar proportions. Many of our sources, however, are fitted with 
significantly lower values between 1 and 4. We interpret these low values to mean that dust in disk 
surface layers has settled vertically towards the midplane. This was expected from the analysis of 
CG97 (see their §3.3), and reinforces similar conclusions by D'Alessio et al. (1999). 

The overall level of infrared excess at A < 100 ^m can also decrease with increasingly edge-on 
viewing angles, as the central portions of the disk are increasingly hidden from view by the flared 
outer "wall" (CG99). Inclination effects appear insufficient to explain the relatively low infrared 
excesses exhibited by sources here, however, because their visual extinctions are < 1 mag and 
their near-to-mid-lR fluxes are not lower than their far-lR fluxes; the reverse would be true for 
significantly inclined sources. Using the results from CG99, we estimate that infrared excesses may 
be depressed by a factor of ~1.5 for our sources due to non-zero inclination; the suppression factor 
due to lower values of H/h is ~2-4. 



4.3. Ice and Silicate Emission Lines 

Observational evidence for silicate emission at 10 ;um from the superheated surface exists for all 
of our sources except AA Tau. In cases where medium-resolution spectra exist (MWC 480, LkCa 
15), the "trapezoidal" shape of the observed emission feature is imperfectly fitted by our model; 
this indicates that actual surface layer silicates have allotropic states (crystalline vs. amorphous) 
and compositions (pyroxene vs. olivine, and Fe:Mg ratios) slightly different from the amorphous 
MgFeSi04 that we employ. Sitko et al. (1999) reproduce the "trapezoidal" shape of the 10 //m 
emission feature using an admixture of amorphous olivine, amorphous enstatite (a pyroxene), and 
crystalline olivine, in roughly equal proportions. The amplitudes of the observed emission bands 
indicate that the dominantly absorbing (emitting) silicates in the surface layer have sizes r < 1 /im, 
well inside the Rayleigh limit. 

Observational evidence for water ice emission at '-^45/um is present in ISO scans of 2 of the 
coolest stars, CQ Tau and AA Tau. The feature at 45 ;um represents the translational mode in water 
ice having the highest oscillator strength. If we examine only the ISO data in the magnified inset 
plots of Figures 9, 10, 11, and 13, and ignore the model fits, there appears to be a trend of increasing 
~45 /im flux relative to ~55 /xm continuum flux with decreasing stellar effective temperature. This 
behavior accords with the trend noted in §3.4 (see Figure 8), whereby the amount of ice present in 
disk surface layers decreases rapidly with increasing T*. 

Notwithstanding this qualitative agreement, the model fits to water ice emission bands (or 
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lack thereof) in the ISO spectra require substantial improvement. The best fitted cases include 
(1) CQ Tau, for which there is even observational evidence of an additional translational band at 
60-65 fiiQ which our model reproduces, and (2) MWC 480, for which the high stellar temperature 
and the relatively small outer truncation radius of the disk suppress water ice emission to levels in 
approximate accord with observations. However, the observed width of the 45 /um band in CQ Tau 
is narrower than what our model predicts. The discordancy of emission line shapes at 45 fim is yet 
stronger in the cases of AA Tau and HD 36112. In the latter case, water ice emission is predicted 
by the model but is not observed. See, however. Figure 14 for an alternative model in which we 
substantially reduce the outer disk radius of HD 36112, thereby reducing the 45 /im flux. 

Accurate reproduction of observed solid-state emission features from water is hampered by 
a number of difficulties. Aside from the crudeness of our 2-layer radiative transfer scheme, these 
obstacles include (1) uncertainties in the photospheric abundances of water relative to silicates (we 
have assumed cosmic abundances with 50% of the oxygen tied up in water and 100% of the iron 
locked in refractory grains); (2) uncertainties in how ice is distributed with particle size (we have 
assumed a constant fractional radial thickness of the ice mantle relative to the radius of the silicate 
core for a power-law distribution of core radii); (3) the probable presence of impurities in water 
ice that can shift band positions and widths; and (4) incompleteness of laboratory data for the 
optical constants of a cosmic mixture of ices in various allotropic states at wavelengths longward 
of 100 //m. Improving the fits by attacking these problems is beyond the scope of our present, 
exploratory work. 

Finally, we note the existence of several apparent emission bands in the ISO spectra that 
we are unable to identify. Most prominent among these is a broad peak near 80 /im in scans 
of AA Tau, MWC 480, and possibly CQ Tau and HD 36112. No resonance at 80;um exists for 
amorphous silicates (Jager et al. 1994); nor is such a resonance measured for the crystalline silicates 
studied by Jager et al. (1998). However, these latter authors also show that peak positions of a 
given vibrational mode shift towards longer wavelengths with increasing iron-to-magnesium content 
(higher effective vibrating masses) . We propose that the 80 fim. feature is caused by a translational 
mode in crystalline olivine having an Fe:Mg ratio intermediate between that of "natural olivine" 
(Mgi,96Feo.o4Si04) and "natural hortonolite" (Mgi.iFeo.9Si04) (Jager et al. 1998, see their Table 
3). If this is the case, we would expect an associated crystalline silicate translational band to appear 
near 51 /um; indeed, such an emission line does appear in ISO scans of AA Tau, MWC 480, and HD 
36112. With regards to these and other perceived emission bands in the ISO data, however, the 
possibility of instrumental error must unfortunately be kept in mind (see §4.2.2 of Creech-Eakman 
et al. 2000). 



4.4. Near-Infrared Excesses 

Our model near-infrared fluxes shortward of lO/um are occasionally underestimated in MWC 
480, LkCa 15, and AA Tau. If we adopt the disk inclination of i ~ 66° for CQ Tau (Natta &: 
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Whitney 2000) and the resultant lowering of its model near-IR fluxes by a factor of 2.5 (see §4), 
then this star's observed near-IR fluxes are also inadequately explained. The relative dearth of 
surface layer emission arises partly because silicate particles are particularly transparent in this 
wavelength regime that exists between the vibrational resonances near 10 /xm and absorption due 
to iron impurities near 1 ^um. It seems possible that thermally-spiked emission from polycyclic 
aromatic hydrocarbons (PAHs) in disk surface layers may help to fill in this transparent region 
of the SED. The strongest resonances are due to C-H and C-C stretching and bending modes 
at 3.3, 6.2, 7.7, 8.6 and 11.3 //m (Draine 1995). Upcoming high spectral resolution observations 
with SOFIA (Stratospheric Observatory For Infrared Astronomy) can test this hypothesis. We 
note, however, that Natta, Prusti, &: Krugel (1993) and Bouwman et al. (2000) argue against 
this explanation for HAe stars AB Aur and HD 163296 based on the need for unrealistically high 
abundances of PAHs and on available ISO Short Wavelength Spectrometer data. 

Other neglected but possibly relevant contributions to near-IR excesses include (1) active 
accretion, which is likely to play an important role at disk radii inside a few AU, (2) reflected 
starlight off the disk surface, and (3) the possibility of dust components in addition to the disk, 
e.g., an optically thin, spherically distributed cloud of silicate/iron grains within a few AU of the 
star (MIVE99; Bouwman et al. 2000). 



5. Summary 

In this work, we have constructed improved versions of 2-layer passive disk models by Chiang &: 
Goldreich (1997). These improvements include an explicit accounting of grain size distributions and 
grain compositions, and numerical solution of the equations of radiative and hydrostatic equilibrium 
under the original 2-layer approximation. We have explored how the SED varies in input parameter 
space and applied our models to observations of 5 T Tauri and Herbig Ae stars. Our principal 
conclusions are as follows: 

1. Hydrostatically flared, passive disks having masses of ~O.Ol-O.lM0 and radii of ~100-250 
AU adequately explain the infrared-to-millimeter wavelength excesses of our sample classi- 
cal T Tauri and HAe stars. Unambiguous determination of the geometry of circumstellar 
dust requires, however, spatially resolved images. Maps from near-infrared to millimeter 
wavelengths generated by the Atacama Large Millimeter Array (ALMA), the Space Infrared 
Telescope Facility (SIRTF), and the Next Generation Space Telescope (NGST) will help to 
break degeneracies inherent in SEDs between, e.g., those of disks and envelopes. 

2. Solid-state spectral features in the mid-infrared (A = 5-60 /xm) appear in emission from face- 
on disks. These emission features arise from "disk atmospheric grains" : grains in disk surface 
layers that are directly irradiated by central starlight. The strongest resonances include the 
10 /im peak from surface silicates at stellocentric distances of a few AU, and the 45;um peak 
from surface water ice at distances of ~100 AU. The strengths of these emission bands relative 
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to that of the adjacent continuum depend on (a) the sizes of atmospheric grains that absorb 
the bulk of the stellar radiation, and (b) the disk viewing geometry. If atmospheric grain sizes 
are within the Rayleigh limit (27rr/A < 1), emission band amplitudes saturate relative to the 
continuum. As atmospheric grain sizes increase beyond the Rayleigh limit (27rr/A -^ oo), 
emission band amplitudes decrease. In addition, as the disk is viewed at increasingly edge-on 
inclinations, emission bands tend to go into absorption (CG99). 

3. Values for Sg, p, rmax,ii Qii and Qq influence the SED most at wavelengths longward of 
100 /um. Their values for a given source, however, cannot be uniquely constrained by the 
millimeter- wave SED alone. One relies on the optical thinness of the disk interior at millimeter 
wavelengths to measure the disk mass in dust: Mdisk oc Fmm oc Tmm oc Skj. The problem is 
that the SED alone cannot disentangle S (which depends on Sq, p, and Oo) from kj (which 
in turn depends on qi and rmax,i)- Spatially resolved millimeter- wave maps can help to break 
the degeneracy between interior grain size and disk mass. 

4. The one disk parameter which appears to be most uniquely constrained by the SED is H/h^ 
the height of the disk photosphere in units of the gas scale height. Its value is roughly 
proportional to the overall level of infrared excess at A < 100 ;um. CQ Tau exhibits H/h = 5.0, 
a value appropriate for gas and dust that are well-mixed in interstellar proportions. Our other 
sources — MWC 480, HD 36112, LkCa 15, and AA Tau — are fitted with significantly lower 
values between 1 and 4. We interpret these low values to mean that atmospheric grains in disk 
surface layers have settled vertically towards the midplane. For standard disk parameters, 
the time required for a 0.1 micron-sized grain to settle from H = Ah to H = Q is S x 10^ yr 
in the absence of vertical gas flow; from H = 4h to H = h, the required time is 8 x 10^ yr 
(CG97; Creech-Eakman et al. 2000). Both these times are of the same order of magnitude as 
the estimated stellar ages. The actual amount of photospheric settling depends also on the 
unknown degree of turbulence and vertical circulation in gas. 

5. Translational lattice modes in water ice appear in emission at 45 ^m and possibly also at 
62 ^m in CQ Tau and AA Tau, two of the coolest stars in our sample (T* < 7200 K). We 
interpret these emission bands as arising from disk atmospheric silicates mantled by water 
ice at stellocentric distances of ~100 AU. The hottest stars in our sample, MWC 480 and 
HD 36112 (T^ > 8400 K), evince no such emission bands. By itself, the dependence on 
stellar temperature of the location of the ice sublimation boundary in the disk surface layer 
is insufficiently steep (approximately asub,s oc T^) to account for the presence and absence, 
respectively, of water ice bands in LWS spectra of CQ Tau and HD 36112; these two stars 
differ in their effective temperatures by only ~15%. 
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A. Solution for Ti{a) and 7(a) 
Equations (4), (5), and (6) combine to yield 
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This is an equation for Ti{a), where 
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Note that (/«j)j is the interior opacity averaged over the Planck function evaluated at Tj; 0, H/h^ 
Tc, T^, and R^ are constants. We rely on the slow and modest variation of 7 with distance to solve 
for Ti(a) in the following manner. We define a logarithmic grid in a = {ai, ... , un} where typically 
N = 300. We begin by guessing a value for 7 at a = ai. This value of 7 is used in equation 
(Al) to solve for both Tj(ai) and Tj(a2) by Brent's root finder (Press et al. 1992). These latter 
values furnish a new 7^7' = 3/2 + (1/2) ln[Tj(a2)/Tj(ai)]/ln(a2/ai) by (A2). This new value of 
7' is then employed in (Al) to compute Tj(a3) and Tj(a4). These, in turn, furnish 7" for Ti{a^) 
and Ti{aQ). Thus the iteration proceeds by updating 7 after every two steps in distance. This 
procedure quickly converges to a smoothly varying solution after the first one or two iterations of 
7. The initial guess of 7 at ai can be improved a posteriori and the calculation repeated. Updating 
7 after every one step in distance introduces numerical instability. That is, if we take 7' to compute 
Tj(a3)|y , then employ Ti{a2)\'y and Tj(a3)|y to calculate 7", and so on, the resultant solution jumps 
erratically with every step. 
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